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Abstract 

I calculate the first corrections to the dynamical pre-exponential factor 
of the bubble nucleation rate for a relativistic first-order phase transition 
in an expanding cosmological background by estimating the effects of the 
Hubble expansion rate on the critical bubbles of hanger's statistical theory 
of metastability. I also comment on possible applications and problems that 
arise when one considers the field theoretical extensions of these results (the 
Coleman-De Luccia and Hawking- Moss instantons and decay rates). 



1 Introduction 



The modern nucleation theory of first-order phase transitions is based on 
the results of Langer in [1] and associated works. These were generahzed 
to quantum field theory [2], at finite temperature [3], and in curved space- 
time [4, 5]. The importance of these considerations is highlighted by the 
fact that the original inflationary models based on a flrst-order cosmological 
phase transition [6, 7, 8, 9] were soon ruled out [10] and replaced by various 
slow-roll models with numerous fine-tuning problems. Of related recent in- 
terest are theories of the landscape and the multiverse [11, 12] for which a 
detailed knowledge of the false vacuum decay rate in curved space-time is of 
vital importance and may lead to a cosmological determination of physical 
parameters such as the cosmological constant. 

Langer's original theory led to the statistical determination of the bubble 
nucleation rate, F, that gives the number of critical sized metastable bubbles 
of the new phase nucleated per unit volume and per unit time, 

where f2 and k are kinematical and dynamical factors respectively. f2 is 
proportional to the physical volume of the system, k is the growth rate of 
the metastable configuration and F its free energy. If is a characteristic 
mass scale of the theory, of the order of the temperature T, then Q and k are 
of order //^ and ji respectively, giving F the correct overall dimensionality. 

In the case of a quantum field theory involving a scalar field <p and a 
Euclidean action functional 



S{ct>) = / 



\{d<t>f + U{<t>) 



(2) 



with a potential that has a relative minimum (false vacuum) and an 

absolute minimum (true vacuum) , the bubble nucleation rate or false vacuum 
decay rate when the field is trapped in the false vacuum is given by [2] 

dN 

Aexp(-S) (3) 



d^x dt 



where B = S^is the Euclidean action of the instanton that is the solution to 
the Euclidean equations of motion and A is a. pre-exponential factor that is 
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given by a functional determinant ratio. If /i is again a characteristic mass 
scale of the theory and A the coupling constant, then Sq is of order 1/A and 
A of order //^. 

If the quantum field theory is considered at finite temperature T that is 
much higher than the inverse critical bubble radius (generically -R* is 

of order l/A/x) then B = S^/T, where 5.3 is the three dimensional action of 
the dimensionally reduced instanton and A is of order [3]. 

The question that I would Hke to address here is what happens to the 
prefactor A when gravitational effects are taken into account, as is the case 
in cosmological applications. It is generally true in flat space-time, and also 
when gravitational effects are weak, that the quantity B gives an exponen- 
tially smaller factor, hence is more important quantitatively. One can, how- 
ever, very well envision a situation, in a landscape or multi verse scenario, 
where both A and B, although still close to their flat space values, have 
an intricate dependence on the model parameters such as the cosmological 
constant, with the net result that the rate T has an actual maximum at the 
observed values. In fact, it is not F per se that is expected to have a peak; 
one is rather interested in suitably defined quantities that measure the rate 
of conversion of physical volume to the new phase [10, 13] like, for example. 



which gives the probability that an arbitrary point in space remains in the 

false vacuum at time t. Here to is the time that signifies the onset of the 
cosmological phase transition, R{t) is the cosmological scale factor, is 
the critical bubble radius at the time of nucleation and V{t) is the velocity 
with which the bubble wall is expanding. In any case F is one of the main 
inputs necessary for the calculation of these quantities. 

The basic problem for the calculation of gravitational effects on A is, 
of course, the lack of a consistent definition of quantum gravity, so, apart 
from dimensional considerations, very few methods have been proposed. In 
[14, 15] the prefactor was calculated for processes involving the creation of 
topological defects, and in [16] corrections to A were estimated with the use 
of the renormalization group. Here I would like to initiate another possibility 
which, although limited in scope from the outset, may be quite useful and is in 
fact, in principle, important, namely the generalization of Langer's original 
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work to the case of curved space-time. Given that the results of [2] and 
[3] are the field theoretical generalizations of the statistical theory of [1], 
it is notable that the corresponding theory, of which the work of [4, 5] is 
supposed to be a generalization, does not exist. The precise formulation 
of such a theory would involve the generalization of thermodynamics and 
kinetic theory in curved space-time and will be the subject of future work. 
Here, instead, I will consider the first gravitational corrections to the pre- 
exponential factor for the previously given description of bubble nucleation 
in first-order phase transitions as given in [17, 18]. That is, I will be working 
in an approximation where the critical bubble radius is much smaller than 
the horizon size and neglecting the gravitational back-reaction, essentially 
assuming that the related Hawking temperature is much smaller than the 
temperature and mass scales of the system. As expected, the leading effect 
of the cosmological expansion will be to increase the bubble nucleation rate 
and an estimate of the corrections is given. One may compare this with works 
that investigate the stability of classical and semiclassical configurations in 
an expanding universe [19, 20, 21] where similar results where found, the 
main difference here being that the bubble configuration considered is already 
metastable. 

Although the result cannot be straightforwardly extended to the case of 
quantum gravity the corrections derived suggest that a more complete treat- 
ment of the nucleation rate is needed; it is important in first-order infiationary 
models, landscape and multiverse scenarios, and may also be of relevance in 
various cases of late time and other first-order cosmological phase transitions 
[22, 23, 24, 25, 26, 27, 28, 29, 30, 31]. 

In Sec. 2, which is essentially the second part of the introduction, I de- 
scribe the results of [4, 5], the problems and the relative corrections expected. 
In Sec. 3 I calculate the first gravitational corrections to the critical bubble 
configurations of Langer's statistical theory of metastability as described in 
[17] and I discuss the results and the approximations involved. In Sec. 4 I 
conclude with some comments. 
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2 Gravitational effects on vacuum decay 



One is interested in a quantum field theory of a scalar field in curved 
space-time with a metric tensor g^i, and Euclidean action 

S = j^g (y^^ + \9,^d'4^d''<t> + (5) 

where TZ is the Ricci curvature, g — det g,j_i,, and for definiteness I will consider 
a potential f/(</') that is everywhere positive and has two minima, as before, 
a relative (false vacuum) and an absolute (true vacuum). 1 will assume again 
that the mass scale is ii and the coupling is A. If the field is trapped in the 
false vacuum and the value of the potential there is £, this will effectively be 
a cosmological constant A = IQnGe and space-time will be de Sitter with a 
Hubble expansion rate H"^ = SirGe/d. From this point on one makes several 
assumptions using the insight from the fiat space-time results described in 
the previous section. First of all one assumes that the solutions to the Eu- 
clidean equations of motion will again describe tunneling and determine the 
exponential factor, B, that is one expects that Euclidean quantum gravity 
and the associated path integral have various features similar to the flat case. 

Then F will determine the nucleation rate of bubbles of the true vacuum, 
which is also de Sitter space-time, with a smaller value of the cosmological 
constant. In terms of the flat case critical bubble radius Rq, which is of order 
1///A, and the horizon radius R^s — ^/H, the results of [4] estimate B in the 
thin- wall approximation as a correction to the flat case value, Bq, 

^ = [1 + {Ro/2Rasr? 
and the radius of the bubble in the presence of gravity, R^, as 

R = -^P (7) 

[l + {Ro/2R,sy] 

where Rq = 3a /e is the bubble radius in the absence of gravity and a is the 
bubble surface tension which is the same for the bubble with gravity in the 
thin wall approximation. 

It is implicit in the CDL formalism that some sort of dilute instanton 
gas approximation must exist in de Sitter space-time in order for the instan- 
ton action to exponentiate, that is one expects an approximation of the sort 
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<< Rds, in which case the gravitational corrections in (6) do not give 
an exponentially small correction to F, in fact they may be comparable to 
gravitational corrections that exist in the pre-exponential factor. The ex- 
pression for B in the opposite limiting case where gravitational effects are 
important, when, for example, the bubble radius is comparable to the horizon 
size, may also be completely different than the Coleman-De Luccia (CDL) 
result, such as is the case in the critical Hawking-Moss (HM) solution [5]. 
Thus the CDL expression has problems in its interpretation in both limiting 
cases. Generally, from a strict point of view, in the weak gravity limit where 
the gravitational corrections in the exponent are of the same order as the 
pre-exponential factor, which has not been calculated, one may say that the 
CDL expression, although natural and reducible to the flat case result, is not 
complete. 

Another curious coincidence arises when one ponders the possibility of 
a thermal interpretation of the gravitational effects, namely when one con- 
siders the fact that de Sitter space-time has a naturally defined Hawking 
temperature, 

T.S - |. (8) 

Provided that a suitable frame is chosen, and a thermal interpretation can 
be given [32, 33], one may expect that, similarly to the transition from A ~ 
Ij,'^ to A in the limit of high temperature, T >> in the fiat 

similar hmit may apply in the gravitational case. This would lead 
to the approximation A ~ i?^ which is expected to hold in other cases in 
the literature [14, 15]. However, the translation of the high temperature 
approximation in the CDL case would read H » A/i, and it is well known 
[34, 35] that when relations like > 4 U" hold, the CDL instanton does 
not even exist. One then may expect that tunnehng is described by the HM 
solution with the pre-exponential factor being approximated by if^, again 
without much justification. 

In view of the possible applications of these results (cosmological phase 
transitions, landscape and multiverse scenarios) one would like to have better 
descriptions and quantitative estimates for F. In principle one expects that 



is the general covariant expression for the production of bubbles of the true 
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vacuum per unit coordinate four- volume, where g = det^f^j,, A has dimen- 
sionahty and is a function of the mass scale, /x, of the theory, the coupling 
constants. A, the Ricci scalar, TZ, and the cosmological constant, A. Also in 
this expression, S is a dimensionless function of the same parameters, and 
the only constraint for both A and B is that they reduce to the flat space- 
time values in the limit of weak gravity. One would expect a relation such 
as (9) to emerge from a saddle point evaluation of a path integral and to 
be meaningful for the cases where the bubble radius is much smaller than 
the horizon; when they are comparable even an interpretation of (9) is not 
straightforward. 

It should be noted that gravitational corrections to the pre-exponential 
factor are expected to exist and can be also estimated in a similar approx- 
imation with an entirely different method that applies the renormalization 
group to the CDL expression [16]. As was mentioned in the introduction, 
here I will try to obtain some more insight into this expression by calculating 
the first gravitational corrections to the dynamical factor k, (of dimension /i) 
that appears in (1), as it was calculated in [17]. The kinematical prefactor fl 
in (1) gives an expression of the form ^f—gy? (times again a dimensionless 
factor that contains gravitational corrections), and similar corrections ap- 
ply in the generalization of the free energy of a metastable system in curved 
space-time. It is assumed that a thermodynamical description of the problem 
can be given in the rest frame of the fluid and the cosmological expansion 
can be treated as a perturbation. This is expected to hold provided that 
a characteristic reaction time, r, for the internal fluid interactions is much 
smaller than . I will also assume that we are in the limits of the thin 
wall approximation in order to treat H as constant in the calculation, and 
also to compare with the CDL results. In principle, however, the equations 
can be solved self-consistently with a variable H. 

In summary, the approximations that will be assumed throughout will be 
that the bubble radius is much smaller than the horizon, that one is within 
the limits of the thin wall approximation, and that relations like r << 
hold. Naturally also the temperature and mass scales of the theory are 
assumed much smaller than the Planck scale. 

The physical situation considered is a metastable relativistic fluid in a 
spatial extend that is large enough in order to feel the cosmological expansion, 
yet small enough in order to treat it as a perturbation. As such, the results are 
not straightforwardly generalized to quantum field theory, they are intended. 
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however, to provide some insight to similar considerations. It should be 
stressed, also, that what is implied is not a breakdown of the semiclassical 
approximations leading to the CDL result, but an interesting dependence 
of the pre-exponential factor on the cosmological expansion rate. This is 
expected for very small values of the cosmological constant A, or R^, << Rds- 
It is suggested that the dependence of F on A (or H) is much more intricate 
than what is described by (6) and the associated expressions of [36]. 

3 First gravitational corrections to Langer's 
theory of metastability 

I will calculate the effects of cosmological expansion described by a metric 

ds^ = -dt'^ + R\t) [dx'^ + dy^ + dz^] (10) 

on the critical bubble solution for a relativistic metastable fluid given in [17] 
which has been used for the description of the QCD phase transition. There 

the dynamical growth rate of the bubble was described by taking into account 
the dissipative effects of a fluid with an energy momentum tensor 

Ti,u ^pgnu+ (p + p)u^ + T^^ (11) 

with 

- (C-2r?/3)(9«ii")(^^. + ii^ii.). (12) 

Here p and p are the fluid energy density and pressure respectively and 
= (1, w) is the four- velocity for a relativistic fluid. As an additional 
phenomenological input one needs the shear and bulk viscocities, r) and C 
respectively. 

The dynamical growth rate, k, is calculated by considering a spherically 
symmetric, exponentially growing, perturbation 

7(f, i) = 7(r) e-^*, (13) 

v{f,t) =v{r)e''\ (14) 



7 



of the also spherically symmetric metastable configuration of a critical bubble 
P = V = When p(r, t) = p{r)+'y{r, t) and the remaining equations are 
used in the equations of motion, self-consistency of the solution determines 

K. 

As explained before, I will calculate the first gravitational corrections by 
determining the effects of the Hubble drag term, H = R/R, on the results of 
[17]. Thus I will assume that a natural free energy description of the problem 
exists in the rest frame of the fluid. Although the considerations here are 
similar to works on stability of classical and solitonic configurations in an 
expanding universe [19, 20, 21], it is important to realize that we are not 
investigating quite the same problem; we are rather interested in the effects 
of cosmological expansion on the actual flat space-time growth rate, kq, of 
the already unstable (metastable) solution. 

The equations of motion 

Vf^T^^'^O (15) 

can be solved approximately in two regions: region (I) which extends from 
just inside the bubble radius i?* to a few correlation lengths, ^, outside the 
bubble surface, and region (II) in distances r greater than the bubble radius 
R* plus a few correlation lengths. The analysis of [17] shows that the fluid 
velocity v{r) has the following behavior: it is very close to zero from the 
center of the bubble up to a few inside the bubble surface, then it rises 
abruptly until in region (I) it starts to fall like 1/r^ and in region (II) it falls 
exponentially to zero. This behavior will be modifled by the Hubble drag 
term, and by matching the two solutions at r ~ R^+ a few ^s, one obtains 
the corrected k, = k,q + 6k,. 

The notation is as follows: overbars will denote the critical bubble solu- 
tion, subscripts will denote the flat space-time values, subscripts I and II 
the respective regions and A will denote the difference of a quantity between 
the equilibrium (subscript t) and the metastable (subscript f) phase. Also, as 
before, ^ will denote the correlation length and R^ the critical bubble radius 
in the presence of gravity. 

One will also make use of the enthalpy density, w — p + p, and the 
bubble surface tension, a, in terms of which the flat case growth rate has 
been calculated as [17] 

M( + 4r^/3) 
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and the kinematical prefactor as 

where the volume of the system, Vol, is usually divided out as in (1). 

One expects these relations to carry over in the presence of gravity, in 
our approximations, by replacing Rq by i?*, keeping the surface tension a 
approximately the same in the thin wall approximation, and the physical 
volume of the system giving a factor of ^^—g. The cosmological expansion, 
however, described by the Hubble drag term, will give an additional contri- 
bution to the dynamical prefactor, that can be calculated from the equations 
of motion. 

The u = equation of motion 

K,^(r) — — - — Ir'^wvir)] — 3Hw (18) 

gives in region (I) 

vi{r) = — + Hr + H^ 19 

and the u — i equation of motion can be simplified in region (II) as 

{k + 2>H)wv{t) = (C + 4r?/3)|; l^^^^yv{T)]^ (20) 

to give 

Mr) e- (21) 



where 



2 3H)w 



The solutions depend on two constants, C and D, the matching is done 
for definiteness at r = i?* + c^, with c a numerical constant of order unity, 
and we get another condition from the fact that, for H = 0, the solution 
should be consistent with the previous result (16). The final result for the 
corrected k — Kq + Sk is 



Sk^KoH^^^^^^[R^ + -]{1 + ^] (23) 
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In order to get an order-of-magnitude estimate for this result one can 
assume that for sufficiently high temperature T one can approximate fjX ^ 
T^/A^, a ~ T^/X to finally obtain the estimate for the corrections to the 
prefactor that arise from the dynamical growth rate: 

(24) 

We see that the leading effect of the cosmological expansion has been, 
indeed, to increase the bubble nucleation rate, and can be significant when 
H ~ A^/^T while, at the same time, H « fi, in accordance with our ap- 
proximations. The surface tension was assumed unchanged in the thin wall 
approximation, the correction to the bubble radius can be estimated from 
(7), it is also subleading, however, in the approximation used. It should be 
noted that this result will also be modified when gravitational corrections to 
Q, the critical bubble radius, and surface tension are incorporated, when one 
goes beyond the limits of our approximations. In any case, even when other 
physical situations are considered, the corrections estimated here also exist 
and can be calculated, for example, by a self-consistent solution of equations 
like (18) and (20). What is more important, and supportive of the arguments 
of the previous section,however, is that the corrections estimated here are dif- 
ferent than what is usually assumed as a naive, dimensional pre-exponential 
factor, for example /x^ or T^, and have an interesting dependence on A (or 
H). 



4 Comments 

The main purpose of this work has been to motivate the suggestion that a 
fuller treatment of the theory of metastability in curved space-time is needed, 
in order to supply the Coleman-Dc Luccia result with possible additional 
gravitational corrections that may provide valuable insight to applications in 
cosmological problems. 

One way to approach this problem is to attempt a generalization of 
Langer's original theory of statistical metastability. The main difficulties 
of this approach stem from the fact that proper definitions of the thermody- 
namical quantities are needed, presumably with the use of relativistic kinetic 
theory in the expanding universe. It was implicitly assumed here that such 
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an extension can be done in the fluid's rest frame and the first corrections 
to the fiat space-time results that were presented show, indeed, the expected 
increase in the bubble nucleation rate due to the cosmological expansion. 
The fuller treatment of the relativistic thermodynamics of first-order phase 
transitions is expected to give additional contributions to the nucleation rate, 
similar to the ones presented here and generally different than the usually 
assumed pre-exponential factor. 
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